Generalized Einstein- Aether theories and the Solar System 
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It has been shown that generahzed Einstein- Aether theories may lead to significant modifications 
to the non-relativistic hmit of the Einstein equations. In this paper we study the effect of a general 
class of such theories on the Solar System. We consider corrections to the gravitational potential in 
negative and positive powers of distance from the source. Using measurements of the perihelion shift 
of Mercury and time delay of radar signals to Cassini, we place constraints on these corrections. We 
find that a subclass of generalized Einstein- Aether theories are compatible with these constraints. 
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I. INTRODUCTION 

The theory of general relativity proposed by Einstein 
explains a wealth of phenomena over a wide range of 
scales. At one extreme, one obtains equations of motion 
for artificial satellites in orbit around the Earth in par- 
ticular and for Solar-System bodies more generally that 
considerably surpass Newtonian gravity in the accuracy 
of their agreement with observations. At the other ex- 
treme we find a comprehensive description of the evolu- 
tion of the Universe and of the growth of structure from a 
spectrum of primordial perturbations which is backed up 
by a handful of precise astronomical observations. This 
apparent success over such a wide range of scales leads 
us to accept general relativity as the correct theory of 
gravitation valid on all scales. 

It is well known that there are a number of imperfec- 
tions in our understanding of the Universe, if we adopt 
Einstein's theory as the theory of gravity and the Stan- 
dard Model of particle physics as the theory of matter. 
For example, there is more gravity in galaxies and clus- 
ters of galaxies than can be accounted for by the matter 
(stars, gas, etc.) that we detect directly through absorp- 
tion or emission of electromagnetic radiation [l|. This 
is explained by invoking the presence of dark matter, 
for which there are several reasonable candidates within 
plausible extensions of the Standard Model of particle 
physics. Furthermore, the Universe seems to be expand- 
ing in a way that general relativity can only accomodate 
by including an exotic form of energy that possesses suf- 
ficient negative pressure to be gravitationally repulsive 
and that currently is the most significant form of en- 
ergy density on cosmo logical scales. This is given the 
name dark energy if it varies in time or space, and the 
cosmological constant otherwise (2]- Finally, in order to 
understand the current homogeneity and isotropy of the 
Universe and to provide an origin for the primordial fluc- 
tuations that grew into cosmic structures, we are led to 
infer the existence of another such form of energy density 
with negative pressure that dominated the evolution of 



the Universe at early times and were the source of the 
quantum fluctuations which were the progenitors of all 
structure. This period of early dark-energy domination 
is called inflation. There are no particularly compelling 
candidates either for dark energy or for inflation within 
extensions to the Standard Model not devised expressly 
for that purpose [3]. 

As remarked, these imperfections can be explained by 
invoking the presence of dark matter and dark energy. It 
is also conceivable that we do not yet have the correct 
theory of gravity and that all or some of these effects are 
due to this fact. This possibility has been the subject of 
intermittent attention for a considerable time @], 
but perhaps especially in the last few years 0, 0]. In 
particular, there has been progress recently in devising 
covariant theories that can, apparently, accommodate all 
these eflfects some level these are, 

consequently, promising alternative theories of gravity. 
As with general relativity, they must therefore satisfy 
the constraints one infers from precise observations of 
the Solar System. 

In this paper we undertake the task of comparing one 
class of modified gravity theories - generalized Einstein- 
Aether theories - to Solar-System observations. In these 
theories, one assumes the existence of a vector field with 
a non-standard kinetic term, and with a time-like vac- 
uum expectation value, at least in the cosmological back- 
ground. Such theories lead to differences from general 
relativistic phenomenology which may be substantial, for 
instance producing modifications to Poisson's equation 
and the Friedmann equation ofprecisely the form posited 
as alternatives to dark matter Q and dark energy [l^ re- 
spectively. They should therefore be constrained by ob- 
servations such as the perihelion precession of Mercury 
and the time delay of radio pulses around the Solar Sys- 
tem. 

The structure of the paper is as follows. We first re- 
view the Einstein- Aether theory in Section |TT1 presenting 
the complete action and the field equations. We then 
focus on the spherically symmetric static metric in Sec- 
tion mil and consider a systematic expansion around the 
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Newtonian solution including terms that decrease more 
quickly with distance as well as terms that grow with dis- 
tance. By considering such a wide range of solutions we 
find constraints on the parameters that define the action 
presented in Section |TT1 In Section [TVl we discuss our re- 
sults. We have organized the text in such a way that the 
main thrust of the calculations are presented in the main 
body of the paper while the details are given in a set of 
appendices at the end. 

II. FIELD EQUATIONS 



Brackets around indices denote symmetrization. 
The equations of motion for the vector field are 

Vo.{:FTp) = 2M^\Aii. (6) 

Variation of the action with respect to A impose on the 
vector field the crucial constraint 

= -1. (7) 

III. SPHERICALLY SYMMETRIC STATIC 
METRIC 



A general action for a vector field A" coupled to gravity 
has the following form: 



S 



R 



IGttG 



N 



+ C{g,A) 



Sm 



(1) 



where g is the metric, R the Ricci scalar and Sm the mat- 
ter action. The Lagrangian of the vector field £{g, A) is 
constructed to be both covariant and local. We use units 
with c = h = 1 and the metric signature is (— , +, +, +). 
Furthermore we demand that the vector field is time-like 
with a fixed length A^Aa = — 1. 

In this paper we will consider a Lagrangian that de- 
pends only on covariant derivatives of A and the time-like 
constraint. It can be written in the form: 



C{A,g) = 



a/3 



cig^'^g^a 



[T{K.) + A(^"v4„ + 1)] (2) 



^Vo.A-'VpA'' 



where Ci are dimensionless constants and M has the 
dimensions of mass. A is a non-dynamical Lagrange- 
multiplicr. 

Note that in the particular case ci — —c^ and C2 = 
we recover the canonical form JC oc Fa/iF"^ where F is 
the field-strength of the four-vector A. 

The gravitational field equations for this theory are 



G 



a(3 



T, 



M'^\A„A 



(3) 



where the stress-energy tensor for the vector field is given 
by 



T, 



1 



a/3 = 7C^a{J^'{J{a Ap) ~ ^'^(q^/3) " J(a/3)^'^)) 



a/3 



-gapM'^T, 



(4) 



where 

ja 



dT 
die 



2/C 



a/3 



,V/3^^ 



ci [{VMiVpA" 



(5) 



(V,v4„)(VM0)] 



In this paper we restrict ourselves to the particular 
case of a spherically symmetric static metric. In isotropic 
coordinates {t,r,6,(j)), it takes the form: 

ds^ = -e^dt^ + e^{dr^ + r^dft^), (8) 

where i^(r) and ^(r) are functions of r only. The vector 
field has only two non-zero components: 



= (A*(r),A'-(r),0,0). 



(9) 



With these restrictions, there are 5 unknown fields, 
v,^,A^,A^ and the Lagrange multiplier A, all of which 
are functions of r only. We need 5 equations to deter- 
mine them. We choose the tt and rr components of the 
Einstein equations ([3]), the t and r components of the 
vector field equation © and the constraint equation ([7]). 
We can then combine these in order to eliminate the La- 
grange multiplier field A, leaving us with 3 dynamical 
equations and the constraint. 

The equations of motion depend on the function T{1C) 
and its first and second derivatives with respect to /C. 
In [13], it was shown that in order to get 'MONDian' 
corrections on galactic scales one has to choose the mass 
parameter M to be of the order M ^ Hq ^ 10~^^ GeV. 
Therefore in the Solar System, where the gravitational 
field is strong with respect to M, the function /C is much 
larger than one, and !F'{1C) can be expanded in inverse 
powers of /C^^^ 



1=1 



JQr/2 



(10) 



where this specific form is suggested by comparison to 
particular modifications of Poisson's equation compatible 
with galaxy phenomenology [12]. We therefore write the 
dynamical equations for a generic term JF'(/C) — 
Note that the leading-order term corresponds to n = 1/2. 
We find 



C" + 2^ + V - — 

r 4 / a2n 



iV(/C)/C (11) 



'jiJC + nhK,') [A-f + (/3/C - nc2e^JC')A^A-' 
+C2e«/C((A'-')' + A'A-"^ + (/4/C + n/5/C') {A'Y 
+ (feIC - ncse^K.') A'A" + c^c^kUa^'Y + A* A*") , 
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\ r 4 2 ) OL2n 2 

+,gi(A^'^)2 + - (ci +C2 + C3)e«(A'^')' 

+ £L_fi(,,')2e-(A*)2 + (ci - c3)i/'e''A*A*' 

+Cle''(A")^ 

and 

= (fti/C + n/t2/C')AM* + (/i3/C + 2nci/C')A*'A'' 

+ (/l4/C-27l(ci+C2+C3)/C')yl^'A* (13) 

+ 2(ci + C2 + C3)/CA'-"v4* - 2ciA:A*"A'- , 
where 

and and hi are functions of ci,C2,C3, the metric 

fields and their first and second derivatives with respect 
to r. The specific expressions for these functions and for 
/C in the metric ([5]) are given in appendix |B] 

Note that equation contains only first deriva- 

tives of the metric and vector fields. It is therefore 
a constraint. For the leading order term at large /C, 
J-'{1C) = the right-hand side vanishes and therefore 
this equation becomes 



A"" = 0, and the motion of the Sun with respect to that 
frame will undoubtedly induce an ^M- Because the 
time-like vacuum expectation value of the aether breaks 
Lorentz invariance, the effect of the solar motion with 
respect to this frame cannot necessarily be accounted for 
by a boost, however it suggests that the contribution to 
wiU be 0{'^sunPsun) ~ 10~^. We defer further con- 
sideration of velocity-induced effects to future work. 

In order to recover Newton's theory we need oi < 0; 
the choice ai = — 1 defines Newton's constant. Post- 
Newtonian corrections are associated with the parame- 
ters 02 and hi. We insert the trial solutions (fT6)l in the 
equations ([7]), PT|) . and and we match the co- 
efficients order by order in — . 

From equation (fTTj) at lowest order, we find: 

nci = 0. (17) 

The case n — Q corresponds to the standard aether th eory 
with T{K.) = /C, and has already been studied [iq^. 
Here we are interested in n ^ which therefore demands 
ci = 0. 

It can be shown that ci — Q leads to no modifica- 
tion to Poisson's equation at all p^. Furthermore, it 
implies the existence of modes which propagate super- 
luminally. There is some debate as to whether this can 
be phenomenologically acceptable [17;]. We therefore try 
here to see whether the inclusion of additional terms in 
the expansion of the metric and vector field can lead to a 
consistent weak field approximation for Lagrangians with 
ci ^0. 



A. Weak field approximation 

In the post-Newtonian parametrization, the metric 
fields are expanded in power of where is the 
Schwarzschild radius of the Sun. This expansion is 
a generalization of the asymptotic behaviour of the 
Schwarzschild metric in isotropic coordinates far from 
the source. Observations of the precession of Mercury 
and the time delay of a signal emitted from the Earth 
and reflected by a satellite or a planet allow us to con- 
strain the first coefficients of the expansion. Hence we 
try the post-Newtonian ansatz for the four fields 

e = 1 + ai h a2 — + ... 

r \ r / 

e« = l + bi^ + b2rA\... (16) 
r \ r / 

r \ r / 

A' = l + ei^+e2f-)' + ... , 
r \ r / 

where ai,bi, di and Ci are free coefficients which are de- 
termined by the equations of motions. 

We shall take A'^ to have no constant component. In 
fact, there is a preferred cosmological frame in which 



B. Additional terms 

In a modified theory of gravity, we expect the metric 
to have additional terms to the ones given in (fTB|) . At the 
scale of a galaxy, the next-to-leading term of the metric 
has to be a positive power or logarithm of r if the mod- 
ified theory is to successfully mimic dark matter on the 
scales of galaxies. More generally, unless there are ad- 
ditional terms that at least decay more slowly than l/r, 
the modified theory will have only minor phenomenolog- 
ical consequences on large scales. Since these additional 
terms have to be completely subdominant in the Solar 
System, their contribution to the equations of motion 
have traditionally been neglected. Nevertheless the form 
of the equations (jlip and (fT2)l suggests that they could 
play a role even in the Solar System. Indeed, in equation 
PT|) the left-hand side is proportional to /C""*"^, whereas 
the right-hand side is proportional to JC. Since /C is ex- 
pected to be much larger than one in the Solar System 
[12 |. the parenthesis on the left-hand side must be sub- 
stantially suppressed. Therefore even if the additional 
terms are small in the Solar System, especially because 
they appear on the left-hand side of equation pT|) or (fT^ 
they may measurably alter the dynamic of the fields in 
the Solar System. Indeed, [H, ^ show that, in modified 
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theories of gravity that seek to explain the observed ac- 
celerated expansion of the Universe, corrections to the ef- 
fective Newtonian potential of a massive body that grow 
with radius are generic. 

Let us take the Minkowski metric as the background 
and add a spherically symmetric perturbation due to the 
Sun, 



= 1 
e« = 1 



-tlj{r). 



(18) 



We now expand the perturbations on a complete set, in- 
cluding both positive and negative powers of r: 



CXD . OO 



(19) 



i=l 



i=0 



Here rgm = (rs/M)^/^. One can show that this is the 
scale at which modifications of gravity occur in theo- 
ries that attempt to modify the Einstein equations to get 



MOND PJ]. For the Sun, r„ 



10^^ km. Since inside 



the Solar System r < 10^ km, we have ^10 ^. 

By a rescaling of the coordinates t and r, we can elim- 
inate Aq and Bq such that the constant term is equal to 
1. The coefficients in front of the additional terms, Ai 
and Bi can be constrained by observations in the Solar 
System. Indeed, observations of the perihelion shift of 
Mercury allow to constrain the parameter 02 to be 0.5 
with an accuracy of Sa2 — 10"'^, while time-delay obser- 
vations allow us to constrain bi to be 1 with an accuracy 
of 6bi ~ lO-'^. 

We derive limits on the coefficients Ai and Bi by re- 
quiring \A,\ (^)' < da2 (^)', where tm ~ 6 x 10^ 
km is the distance between Mercury and the Sun. The 
same argument applies for the other metric field. We find 

l^^il (7^) ~ ^bi^. The best constraint on hi comes 

from the Cassini satellite when it was at ~ 10^ km from 
the Sun [23. 

These limits would be exact if the only non-zero term 
were the one under consideration. Otherwise, the true 
limits may be either weaker or stronger (as the sum of 
a series can easily be smaller than some of its terms), 
or even have the opposite sense (> instead of < or 
vice versa). For example, a correction of the form 
exp(— Cr/rgm) with C > requires a lower limit on C 
rather than the upper limits that would be derived from a 
term- by-term analysis. Thus the limits we shall proceed 
to derive from a term- by-term approach will be sufhcient, 
but not necessary. 

Indeed, we show that a set of consistency conditions 
between the coefficients of the powers of r/rgm can be 
satisfied simultaneously. Thus an acceptable solution ex- 
ists which has an expansion that includes both negative 







A\A'' 


Ao = 


Bo = 


|i?o|,|i5o| <3-10-» 


l^il ^ 10-"* 


\Bi \ < 10-" 


\Di\AE^\ < 10-« 


\A2\ < 10-" 


|B2|<5-10-« 


|752|,|£2|<5-10-4 



TABLE I: Constraints from observations on the coefficients 
of positive powers of -;r—- 



and positive powers of r/rgm- Future detailed study of 
the Schwarzschild metric in Einstein- Aether theories may 
yet identify an expansion basis that is better adapted to 
the physics than polynomials in r. 

One might also ask if these extra-terms in the metric 
are physically acceptable if considered in isolation. After 
all, they diverge when r goes to infinity. However, the ex- 
pansion pU]) is itself appropriate only for r <C Af-^, and 
in this paper we are looking to find solutions consistent 
with observational constraints in this range. Meanwhile 
spatial infinity lies in the region r 3> M~^, where an- 
other limiting form for T is appropriate. We discuss the 
expected geometry in this region in Appendix VK\ 

Through the equations of motion, the perturbations of 
the metric created by the Sun lead to perturbations of 
the vector field: 



A''{r) = a{r) 
A\r) = l + /3(r). 



(20) 



We can again expand these perturbations on the com- 
plete set: 



(21) 



i=l 1=0 
00 - 00 / \ » 

2—1 i—O \ y / 



In the following we restrict ourselves to the leading 
order term n = 1/2 in the expansion (jlOp for J-' . In ap- 
pendix [C] we show that a(r) and l3(r) are of the same 
order-of-magnitude as the perturbations of the metric 
(j){r) and ip^r). Since the leading term of the metric in the 
Solar System is ^, we will take this value at the edge of 
the Solar System as a limit for the peculiar terms. This 
means that 



\D, 



gin 



< 



and \Ei 



gm 



< 



(22) 



In table m we summarize the resulting constraints on the 
coefficients Ai,Bi, Di and Ei with i < 2. 

Next we wish to use our expansion of the metric pop 
and the vector field (|2ip to solve equations (jlip . (|12p . 
P3p and the constraint equation ([7]) order- by-order in r. 
We first rewrite the positive powers of r in the form 



Ar. 



gm 



A.. 



gm 



Ar. 



(23) 



5 



where An = An (7^) ■ We do the same for the co- 
efficients of the other fields: i3„,i?„ and Since 
~ 10""'^^, the constraints on the An in table |T] be- 

come much more stringent constraints on the An (see 
table mi in appendix ID 2p . 

We use a perturbative expansion in the coefficients. 
Moreover we assume the hierarchy 1 ^ \Ai \ :$> IA2I ^ ... 
for each of the four fields in order to truncate the number 
of terms present at each order. The detailed calculation is 
presented in appendix |D] Here we summarize the results. 

Equation leads to two possibilities. The first is 
that 4 = and ^ for aU s, i.e. A'ir) = 0. 
This means that the only degree of freedom of the vector 
field is ^* (r) which is then completely determined by the 
metric through the unit time- like constraint. We do not 
study this specific case in this paper, but rather concen- 
trate on the other possibility A'^{r) ^ 0. Nevertheless 
the former solution could be relevant in some particular 
cases. For example in [2l|, it has been shown that in the 
usual Einstein- Aether theory, corresponding to n = 0, 
there are regular perfect fluid stars with a static (i.e. 
A'^{r) = 0) aether exterior. 

The choice A'^{r) ^ implies that the parameters ci , C2 
and C3 satisfy the constraint 

(C1+C2 + C3)S^-3(C1+C2+C3)S-2(C1-C2+C3) = , (24) 

for some positive integer s which is such that di — for 
all i < s. 

We have explicitly calculated the case s = 1, which 
corresponds to the constraint C3 = — ci. From the posi- 
tive orders in — we find that ai is unconstrained. If we 

r ^ 

set ai = — 1 to recover Newton's theory, we find 



61 = 1 + 0(10"^) and 



02 = ^ + 0(10-9) , 



(25) 



which is in complete agreement with the observations 
which measure 



bf" = 1 ± 10~ 



and 



^obs 



(26) 



Moreover, we find that ci no longer has to be equal to 
zero. We have four solutions for ci and C2 which are 
negative and therefore causal p^ . They are given in 
appendix lD 41 In each solution, ci and C2 are of the order 

of ( J^) ^. Since ( J2_) ^ 10-18, either ai must be 

very small, or ci and C2 are very small. 

In figure [1] we plot gtt{r) = —e'^^^^ and grr{r) = e^^'') 
(solid line) for the values of ai to a4, bi to 64, Ai, A2 
and Bi calculated in appendix ID4bl We find that B2 
is unconstrained by the equation of motion, except that 
IB2I < 5 • 10^9 (see table |l|. We choose to saturate this 
maximum value for the plot. We also plot gu and g^r 
from General Relativity, i.e. the Schwarzschild metric 
up to the order (^) (dotted line). We see that inside 
the Solar System (r < 10^ km) the two models are equiv- 
alent, but at a scale rm — 10^^ km they start to diverge. 
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FIG. 1: The metric fields gu and g^r as functions of r, in 
the aether theory (solid line) and in general relativity (dotted 
line). We see that at the scale r™, ~ 10^^ km the two models 
start to diverge. rg,„ ~ 10^^ km is the scale at which all 
growing corrections are expected to start to dominate. 



The Schwarzschild metric tends to gtt = and grr = 1, 
whereas the aether metric starts to grow. The scale r„i is 
the one at which the first two growing corrections (r and 
r^) become dominant. The subsequent growing terms 
could start to grow earlier. Nevertheless the limits set 
by hand on the coefficients Ai and Bi (see table |T| en- 
sure that they cannot become important inside the Solar 
System. Moreover — 10^^ km is the scale at which 
modifications of gravity should occur. Hence we expect 
the growing terms to remain small for r < rg,„. Nev- 
ertheless, only the full calculation of all the coefficients 
could confirm this limit. Note that gtt{f) passes through 
zero at r ~ 10^^ km and is therefore singular. However 
the series cannot be taken seriously for r > rgm, 
where they may not converge. 

Finally we have determined the structure of the equa- 
tions at each order. For negative powers in r (s > 1), 



S = 



s 



( -1 





( 




( ^ 


bs 










-"3 









where 



(27) 



dis{s—l) 





s{s-l) 
— s 




-2 \ 





ci(s-l)(s-2) dis{s-l) 
4 2 



(28) 



(s) 

and HI"' are functions of the previous coefficients - 
ai, ...as_i,6i, ...,6s_i,di, ..,ds-i,ei, ..,es_i. Therefore 
at order s a unique solution (a^., bg, d^, Cg) exists if 



detS": 



-ys'(s-l)'(s-2)^0. (29) 
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Since ci ^ 0, this condition is satisfied for s > 2. The 
order s = 2, for which the determinant vanishes, has been 
solved explicitly. The solutions are not unique and are 
given in appendix ID 41 For each order s > 2, we can find 
the unique solution as a function of the lower orders. 

The same structure repeats for positive powers of r 
(s > 2): 







( H^^ 




Bs 








Ds 
















1 











-s(s + l) 






cis(s+l)c/i 



V 





[c2{U-s)dl- 
ci(s-l)(s-2)i 



where 



2 ^ 





c.is{s+l)di 



(30) 



(31) 



(s) 

and HI are functions of ai, 6i, di and ei, and 
of the previous coefficients - Ai, Ag-i, Bi, Ss-i, 
Di, Ds-i,Ei, ...,Es-i- In appendix ID 21 we argue that 
these functions are at most of the same order of mag- 

" (s) 

nitude as the coefficient of order s. The Q] are func- 
tions of the coefficients As+i, As+2i -Bs+i, Ds+i, 
Eg^i, .... Generically the are much larger than the 

^ ( s) ~ 

Q\ since the constraints on the coefficients A^ are much 
more stringent for larger s. We can, in general, therefore 
neglect the Q^'^ (See appendix ID 21 for a more detailed 
discussion) . 

A unique solution (A^, B^, Ds,Es) exists then at order 
s if 



dots' = -2s^{s- 



)(c2(14- 



.s)dl-^{s-l){s-2] 



^0. 

(32) 

The order s = 0, 1 and 2 have been calculated in appendix 
lD4l For each of the solutions (|D42l) to ([D45)) we find that 
det S* ^ V s. 

The structure of equations (^7)) and ([50]) shows that 
there exist solutions to the equations of motion that can 
be expressed as the expansions ([TS]) for the metric and 
(|2T|) for the vector field. Indeed, since at each order the 
determinants det S and det S are non zero, we are en- 
sured that a solution exists. We conclude that it is not 
sufficient to calculate the first two orders, even if we are 
only interested in the values of the Post-Newtonian pa- 
rameters. In order to be sure that we have a bona fide 
solution of the equations of motion and that we will not 
encounter an inconsistency at any given order, we need 
to calculate the determinant of the system at every order. 

IV. CONCLUSION 

In this paper, we have studied the constraints we im- 
posed from Solar System observations on generalized 



Einstein- Aether theory. We have considered an expan- 
sion of the metric around the Newtonian solution includ- 
ing the usual Post-Newtonian terms, which are negative 
powers of the distance r from the Sun, but also terms 
increasing with the distance (positive powers of r). The 
aim of this complete expansion is to take into account 
in the Solar System the effect of modification of gravity 
at large (Galactic and cosmological) scales. These effects 
are usually neglected in the Solar System and one con- 
siders only the decreasing terms in the metric expansion 
^15, 16, 23]. Nevertheless, as long as the increasing terms 
are sufficiently small to be consistent with observations, 
nothing forces them to be completely absent in the Solar 
System. Moreover in the case of generalized Einstein- 
Aether theory, we found that the increasing terms play a 
crucial role, since without them the theory suffers from 
acausality. Indeed, if we neglect them, the equations of 
motion imply that one parameter of the theory ci — 0, 
leading to superluminal propagation of the aether field 
perturbations. If we consider the full expansion, the 
pathology disappears. 

Moreover, we found that the Post-Newtonian parame- 
ters related to the precession of the perihelion of Mercury 
and the time delay of radio pulses are in agreement with 
observations. We found also constraints on the increasing 
terms coming from these observations, and constraints 
on the parameters ci, C2, C3 and ai from the equations of 
motion. Indeed, whereas the ai represent behaviour in 
a particular regime of the theory (i.e. quasistatic config- 
urations where the gravitational field is typically much 
larger than the mass scale M), the Ci are parameters 
of the Lagrangian which affect all solutions. It may be 
shown |12l . |2J| that simultaneous consideration of very 
weak gravitational fields, the requirement for a realistic 
background cosmology, and suitable growth of large scale 
structure will favour the Ci to be 0(1). We may imme- 
diately see then from equations (|D42p - (jD45P that this 
corresponds to an upper limit on the number ai of 10^^ 
to 10^^. It is interesting to compare this bound to that 
obtained by considering the influence of such a term on 
the Poisson equation then deducing the perihelion pre- 
cession due to the extra force provided by the modified 
potential in the context of Newtonian gravity. This was 
done in [4i| where it was found that the resulting perihe- 
lion shift S(j)a of Mercury's orbit per revolution, to lowest 
order in eccentricity, was given by: 



lOTrMai — 



(33) 



where r is the semi-major axis length of Mercury's or- 
bit. Taking M = 1.2 x 10~^cm/ , this yields a preces- 
sion per revolution which is approximately ai/5 of the 
prediction due to general relativity. Therefore, as ob- 
servation agrees with the prediction of general relativity 
to within ~ 10~^, it was argued, ai must be less than 
around 5 x 10~^, therefore several orders of magnitude 
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greater than our analysis will allow. Clearly then, a ben- 
efit of recovering the modified Poisson equation from a 
set of generally covariant field equations has been to al- 
low a broadening of the scope of analysis of the theory's 
implications in the solar system and in this case, in con- 
cert with additional constraints placed on the theory's 
parameters in other regimes, this has allowed for a sig- 
nificantly more severe restrictions on the permissable size 
of such a modification. 

Furthermore, we have developed a general method to 
test that the metric expansion is a solution of the equa- 
tions of motion. The constraints that we obtain must be 
interpreted, in part, as consistency conditions on the the- 
ory. Indeed, it is not sufficient to calculate the first coeffi- 
cients in the expansion to confront the theory with Solar 
System observations, but one has to study carefully the 
structure of the equations at each order to be sure that 
no inconsistency will invalidate the results. Our results 
show that, even though we have an infinite hierarchy of 
equations, these are solvable in terms of the fundamental 
parameters of the theory. Hence, generalized Einstein- 
Aether theories are viable theories of gravity within the 
Solar System. 
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APPENDIX A: ON ASYMPTOTIC FLATNESS 

In the range ?' ^ we expect the geometry to pro- 
duce MOND behaviour which we will define as follows: 
We can decompose the 4-velocity of a test particle 
into a part proportional to the velocity of a static 
observer in the spacetime and a component satisfying 
C'^Vf_i = 0. Normalization enforces U^^ = {l+v'^)^ C^+yf^ 
where v'^ = g^^v^v'^ . The MOND regime is the regime 
for which v^{r) = (rsM/2)2 . In the weak field limit the 
constancy of at large radii can be interpreted as the 
levelling off of orbital velocity for objects in circular or- 
bits, for instance the so-called 'flat rotation curves' of 
stars in approximately circular orbits far from the centre 
of galaxies. 



We would like then to obtain an insight into the ex- 
pected geometry the model describes in this regime. 
It can be shown that another relativistic theory which 
allows MOND, Bekenstein's TeVeS theory, can be re- 
written in a form very similar to the model discussed 
here (25| i.e. it may be written as a theory with a single 
metric accompanied by a derivative coupled Lorentz vi- 
olating vector field with noncanonical kinetic terms. In 
TeVeS for the static spherically symmetric case where the 
vector field points in the time direction the metric for 
r ^ takes the form in Schwarzschild coordinates 

In 

ds'^-i ^ + + p'dVL^ , (Al) 

VPo/ {\-nY 

where po is a constant of integration, n = psMj (4 -|- 
PsM) <C 1 and ps is the Schwarzschild radius of the Sun 
in Schwarzschild coordinates. The form of goo is fixed, 
via the geodesic equations, by the requirement that is 
independent of p. It can be shown that this metric has 
nonvanishing components for the Riemann tensor i?"^^!*- 
= [n - l)ns\v?e,R'''g^g = -n(n - 2), i?*^,^ = 
—n{n — 1). It is this tensor, for instance, which is a mea- 
sure of the change AiS^ in components of an arbitrary 
one-form Sa parallel transported around a small closed 
curve at some point p in the spacetime [20| . Explicitly: 

^Sb^^R\,^Saj>x''dx^ (A2) 

Clearly this is not generally zero if, as in this case, 
^""bcd does not vanish and so we may conclude that the 
spacetime described by the above metric is not flat. Due 
to the similarity between the aether stress energy tensors 
in the two models, we expect the metric to take a very 
similar asymptotic form for the model considered here 
when the aether points in the time direction. However, 
though part of the geometry (and thus part of the neces- 
sary behaviour of the aether) is fixed by the requirement 
that the MOND regime exists, a more thorough treat- 
ment is required to see how this is affected by allowing 
for a radial component to the aether. It is interesting to 
note that for the case T — JC, asymptotically flat solu- 
tions were only found for cases where the aether aligned 
with dt .21]. It is tempting then to conclude that for mod- 
els with forms of T allowing a MOND regime for A cx 9t , 
the addition of radial terms, if allowing a MOND regime 
at all, would not reduce the asymptotic curvature and 
that the absence of asymptotic flatness is generic. 

It may be checked that when (|Aip is written in 
isotropic coordinates, the metric component goo varies 
with the isotropic radial distance r as — r^" whilst grr 
varies as r~^". The approach to the weak field limit 
of the metric can be found by requiring that goo is 
rather close to one, so we may expand — (r/ro)^" = 
-exp(2nln(r/ro)) ~ (H- 2 ln(r/ro)) = -(1 + 2$). The 
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difference in the Newtonian potential <I> from the onset 
of MOND at the typical scale of a galaxy to the current 
Hubble radius is of the order 10~^ and so though the 
metric components do asymptotically diverge, they are 
well within the weak field regime within the current cos- 
mological horizon. It is interesting then that the metric 
is approximately Minkowskian for r ~ M^^ but diverges 
to greater curvature at far shorter and far greater dis- 
tances from the source. However, the mildness of the 
deviation from flatness even on cosmological scales ren- 
ders it unclear whether the absence of asymptotic flatness 
of the kind in (jAl|) would have observable consequences. 
It is also unclear whether models producing MOND and 
a vector aligned with dt are more stable than solutions 
admitting a nonzero radial component to the aether or 
whether they do tend to arise preferentially from astro- 
physical initial conditions. Indeed it has been claimed 
that the static spherically symmetric solution with zero 
aether radial component in TeVeS is unstable in the solar 
system with respect to certain perturbations, the insta- 
bility, partially manifest as an increasing radial compo- 
nent of the aether, growing on unacceptably short time 
scales [13 . We postpone more detailed analysis to future 
work. 



+C2 
92 = 



C3 /4 



(Cl+C3)(^ + — + - + ^ 



8 9^'^ 4(1/' + 3^') 



/ 8 V"- 94'" 



4 



(ci+C3)C' + C2(i/' + 3e' + -) 
4 2(z/' + e') 



/ii = -cAv''' 

\ r 

/ , 4 2{v' -i!) 

-t-ca \v i 2 ^ 

\ r r 



v" - C 



II I o fll 



^ -ci (Si^' + + -) and 
hi = (ci+C2 + C3)(i/' + 3e' + ^) 



APPENDIX B: SPECIFIC EXPRESSION OF THE 
EQUATIONS OF MOTION 



APPENDIX C: ORDER OF MAGNITUDE OF 
THE VECTOR FIELD PERTURBATIONS 



The explicit expression for /C in the spherically sym- 
metric and static metric is 



C2 



r r 2 4 4 

4 



(ci + C3)C' -I- C2 (- + I/' + i^'y^A-A-' (Bl) 



,12 



+ e''-^{c:i~c^)[—{A'f + v'A'A' 



+ (ci+C2 + C3)(A'-')2-cie''-«(A*') 

The functions appearing in the equations of motions 
(fTT|l . (fT2l) and ^ are given by: 



(Cl +C3)( 



3j/'e' 2u' 



(B2) 



+^^[-^ + + 2 



h - 



2 \ r 



h = y (^3 -cO(— + — + - + -' 



In this appendix, we compare orders of magnitude of 
the different field perturbations. We know that in the 
Solar System the metric perturbations are of the order 
(j){r) ~ ■0(r) ~ — . We want to calculate the order of mag- 
nitude of the vector field perturbations a{r) and /3(r). 

At first order in the perturbations: 



(Cl) 



A'' = 13' and A''' = a' . 
The constraint ([7]) gives 

-(l + 0)(H-/3)2 + (l + Vj)a2 = -1 . (C2) 
At the lowest order in each perturbation we have 

-<j)-2l3 + a^ = . (C3) 
We have four possibilities: 

1. The three terms are or the same order of magnitude 
(3 ^ ^ 4) which implies — 4) + 2(3 \ 

2. <^ (f> ^ (3 which implies ~ —2(3] 

3. (3 <^ (j) ^ which implies ~ (/); 

4. (p ^ (3 ^ which implies ~ 2(3 . 
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We see directly that case |3] is not possible, since <f) = 
—Tsjr < at first order in the Solar system. We use 
equation (fT3|) to exclude the possibilities [T] and [H In- 
deed, in these two cases we have a ^ f3,4>. Using these 
constraints we find 



1 

M2 



2a (ci + 2c2 + C3) 4a'ac2 



(C4) 



+ (Ci + C2 + C3)a'2 



and equation becomes 



= 2(ci+C2+C3) 



2^ + 2— + a'' 



/C 



C2 h (Cl + C2 + C3)a' 

r 



/C' = 



(C5) 



The only solution of this equation is a(r) — 0, which 
implies that equation (fT3|l doesn't allow the spatial com- 
ponent to be much larger than the other perturba- 
tions. So the only possibility is that the dominant con- 
tribution to the time component of the vector field is 
/9 = and that the dominant contribution to the 

space component satisfies a <C — . So it is legitimate to 
assume that also a '--^ — at first order. 



APPENDIX D: RESOLUTION 

In this appendix we present the detailed derivation of 
the solutions. The aim is to take the expansion p9)l for 
the metric and (pij) for the vector field, to insert it in 
the equations of motions and to solve order by order in 
power of r, using a perturbative expansion for the posi- 
tive powers of r. The orders r^^ and r^^ give the values 
of the post-Newtonian parameters bi and 02- The or- 
ders r~^,r~'^ and r°, and allow us to constrain 
the values of the Ci which are compatible with the data. 
Nevertheless it is not sufficient to calculate only these 
orders. Indeed, each equation generates an infinite num- 
bers of positive and negative orders, and the expansions 
(fro)) and (PT|) are solutions of the equations of motion 
only if each of these equations has a solution. Therefore, 
it is crucial to understand the structure of the equations 
at each order and to verify that one introduces a suffi- 
cient number of new coefficients at each order so that 
the equations do not lead to constraints between coeffi- 
cients which have already been determined, leading to a 
possible inconsistency. 

In ID li we consider only the terms with negative powers 
of r, that means the terms which appear in the usual 
Post-Newtonian parametrization . In ID 2l we will consider 
the positive powers. Finally, in ID 31 and ID 41 we apply our 
result to equations mi), ([121), and 



1. Negative powers of r 

We have four fields A'',A*,e^ and e'^, satisfying four 
equations. Three fields, A*,e^ and e'^ have an expansion 
of the form 



(Dl) 



where x = The other field has no constant term in 
the expansion 



(D2) 



In the following we restrict ourselves to the two fields: 
and e"^ and two equations. The generalization for the 
four fields will then be straightforward. 

We assume that the two equations governing the two 
fields have the form 

m 

J2f^r'"'ieT'{e'"f'{e'"'^^{A'^ f^{A'^'Y'{A"')P' =0 , 
1=1 

(D3) 

where a' stands for ^ and the exponents are all positive 
integers or null, except ct^ which is negative. Actually 
ai will in principle be negative, since ly' = {e'^)' {e'^)^^ 
and ly" = {e")" {e")-^ - {e^' f {e")-^ , but we can always 
multiply by (e")-' so that is a positive integer or null 
for all i. 

The form (|D3p is not exactly the one which appears in 
equations (fTT|) and (fT^ . because of the term K^, where 
n is an integer or a half integer. Nevertheless, we will see 
later that we can rewrite the problem in such a way that 
the equation is the form (|D3|) . 

Since all the terms in the sum must have the same 
dimension, we have that —at + (3i + 2"fi + Vi + 2pi = c is 
the same for all i. Therefore, we can solve for ai, and then 
multiply every term by r'^. We insert then the expansion 
(IDll) and (ID2D. and their derivatives to obtain 



J2 /.(-1)'3-+''-x-('3.+27>+-.+2p0 1 + ^ a„a;" 

i—l \ n—1 / 

Coo \ / 00 \ T» 

/ 00 \ A". / 00 \ 



/ 00 



n{n + l)c?„x" 



+2 



(D4) 



Each term of the product can be expanded in power of x 
using the binomial theorem 



s1"= i + E«"^" =1 



(D5) 
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+ + • • • , ds-i)]x^''+'-^ + ... 

where gs{di, • • • , dg-i) is a function of the coefhcients of 
order 1 to s — 1 . 




(D7) 

and the same for A^" . 

We can insert these developments in equation (|D4p and 
solve order by order in x. The lowest order is 

(D8) 

where Oi = |2 G {1, m} such that /9i + 7i + /ii + f ^ + 

Pi = (5 is the smallest exponent of the sum|. 

The next order in x contains two different contribu- 
tions. The first contribution comes from the terms in 

01 where we take the next order in the expansion for 
one term of the product, and the lowest order for the 
others. This contribution contains two new coefficients 

02 and d2 which appear linearly, and also the previous 
coefficients ai and di. The second contribution comes 
from other terms in the sum, i G O2, where O2 = {* ^ 

{1, m} such that /3j + 7^ + /ij + t',; + = 5 + l|. This 
contribution contains only ai and c?i. The next order 
x^^^ is therefore 

/.(-l)'^-+'^*2^-+''-{(2A + 37.)«f +^'-'<''"^''''"a2 
+ + 21,, + 3p,)af +'^'+'''-id2}x^+' 



ISO 



+H\^\ai,di)x 







(D9) 



More generally, at the order 2:''+'' ^ (s > 1) the two equa- 
tions have the form 

F['\ai,di)-as + G['\audi)-ds 
= H[''\ai, ■ ■ ■ ,as-i,di, ■ ■ -,4-1) , 

F^'\ai,di)-a, + Gi'\ai,di)-ds 

= i/^'^(ai,---,a,_i,di,---,d,_i), (DIG) 

where 

F^'\ai,d,)^ Y /,(-l)^'-+''-2^'+''--i. 
leOi 

X (2sA + sis + l)7.)af+^'''rff 



leOi 

X (2^, + 2siy, + s{s + l)p,)af +"•+''■-1 , 

(Dll) 

and _ffp^(ai, • • • , Os^i, di, • • • , ds-i) is a function of the 
coefficients of order less than s. F2*''(oi,di) and 
G2^\ai,di) have the same form but with different fi and 
Oi- 

A unique solution (os, ds) exists at order s if 

,^jFi^yd,) G^\a„d,)\ 

\F^'\a,,d^) Gi'\a,,d,)) )^ ^ ' 



Since the functions -Pj'*'' and G^*'' contain only oi, di 
and the parameters /; for i G Oi, it is easy to calculate 
them at each order, without solving the entire system. 
We only have to solve explicitly the first order. Then 
if the determinant is non zero for all s > 1 we know 
that a unique solution exists which can be determined 
as a function of the previous coefficients. On the other 
hand if the determinant vanishes at some order s, an 
inconsistency may occur. Indeed, equations (|D10[) can 
become a constraint between the previous coefficients and 
lead to a contradiction. Therefore, if one wants to test 
the validity of a theory in the PPN parametrization, it 
is not sufficient to calculate the first coefficients in the 
expansion and to compare them with the observations. 
One has also to calculate the determinant ()D12|) at each 
order s, to insure that no additional constraints on the 
first coefficients will occur from higher orders. 

The generalization of this method to four equations 
with four fields is straightforward. Each field with no 
constant component generates a function of the form 

( s) 

, whereas fields with a constant component gener- 



ate a function of the form G 



2. Positive powers of r 

In this section, we add positive powers of r to the ex- 
pansions (|DT]1 and (|D2|l : 



n—l ^ ^ ^ n— 1 

00 00 



n=l 



n=l 



(D13) 



where e = ~ 10 . Equivalently for 

Tgrn. 

OO - oo 



n=0 
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These new terms add an infinite number of contribu- 
tions to each previous order x^~^^ . Furthermore they gen- 
erate new lower orders x^~^. To simpUfy the problem 
we must take into account the fact that the coefficients 
in front of the negative orders A„e" (respectively £>„e") 
have to be small in order not to be observed in the Solar 
System. Indeed the constraints of table H] are equivalent 

to |^„|e" ^ da2 (7^7) ^1- Therefore, we use a per- 
turbative expansion for A„e" and DnC^- Furthermore 
since ~ 5 • 10^^ ^ 1 we have the following hierarchy: 

1 > |Ai|e > |A2|e^.... (And similarly for the So 
we start studying the effect of Dq (remember that Aq has 
been put to zero by a rescaling of t and r), then and 
and so on. And we neglect all these new terms with 
respect to 1. In other words, we take them into account 
only when they introduce a new order to the equation 

Let us define new coefficients to simplify the formulae. 

i„ = e"A„ £)„ = e"i^„ . (D15) 

In terms of these coefficients, the constraints coming 
from observations are as follows: 





e« 


A\A'' 


io = 


Bo = 


|A)|, l^ol <3-10-3 


|ii| < 10-2^ 


\Bi \ < 10-22 


pil, iBil < 10-" 


|i2|< 10-^2 


\B2\<5-W''' 


1^21, \E2\ < -W'^" 



TABLE II: Constraints from observations on the new coeffi- 
cients . 



a. Effect of Do 



(D16) 



The term Dq introduces a new order x 



s-i 



b. Effect of Aix ^ and Dix ^ 
At first order in Ai and Di we have 



5-1 



(D17) 



1 



X 



ji=i 



ai h 1 -|- aittix ■ 

X 



(D18) 



(A^ 




(D19) 



13. 



+ (-l)^-iA(A-l)2iiaf-2a2x2'^--i 
-f(-l)'^'/3.afx2^-+-.. , 



(D20) 



and cquivalcntly for A^' . e"" and A'^" remain the same 
as previously. Therefore, we see that Aix~^ and Dix~^ 
introduce a new order x^~'^ and also a contribution to 
the order x^~^ . 



c. General: AsX " and DgX 



More generally, the term AgX " introduces a low- 



est order 



„<5-s-l 



linear in As- The following terms 



As+\, As+2--- also contribute to the order x^~'^~^ , but 
there are negligible with respect to Ag. Furthermore, the 
previous terms As-i^ As-2--- also appear at order 
but non-linearly. We can show from the observational 
constraints that we have on the coefficients, that these 
non-linear contributions can be at most of the same or- 
der of magnitude as As, but not larger, since they contain 
a term of the order e'* . Of course the same kind of terms 
are introduced by DgX^" . 

Therefore, the equations at order x''-'*-^ (s>0) are 

Fi'\ai,di)-As + G^{\ai,di)-Ds 



H{'>{Ai,...,As-i,Di,...,Ds-i,ai,di) 



F^'>{ai,di)-As + G'^>{ai,di)-D 

+0^^^(i,+i,...,^,+i,...) 



(D21) 



where 



Fr{a^,d,)^Y. /.(-l)^-+'''-i2^-+'''-i 
(2sA-s(s-l)7,)af+^'-idf+'^' 



+Pi 



X 

G['\a,,d,) ^ /,(-l)ft+''-i2^-+'"-i 
leOi 



X ( - 2/i, + 2sv, - s{s - l)p,)af +"'+"'-1 . 

(D22) 
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ijj^-* can be at most of the order of magnitude of Ag 
and Ds whereas Qj are much smaller since they are of 

the order of magnitude of ^s+i and Ds+i- 

A unique solution {Ag, Ds) exists at order s if 



det 



7^0 



(D23) 



In this case, we can determine the solution at order s as 
a function of the previous coefficients and the following 

" (s) 

coefficients. Generically the contribution from the 
previous coefficients is much larger than the one from 
the following coefficients Q\^\ since the constraints on 
the previous coefficients is less stringent (see table |lTl . 
We can therefore neglect the following coefhcients and at 
each order determine the solution as function of the pre- 
vious coefficients. Nevertheless for some specific case the 
contribution from the previous coefficients can become 
very small and even vanish. In this case, we have to take 
into account the contribution of the following coefficients. 
This means that the solution {Ag, Ds) will be of the order 
of magnitude of {As+i, Ds+i). At order s-l- 1 we can then 
neglect As and Ds in Hj^^~^^^ which appear non-linearly 
and are therefore much smaller than the order s + 1 lin- 
early. This means that we can determine (As+i, Ds+i) 
as a function of the previous coefficients, but neglecting 
the order s. If the determinant (jD23p is non zero for all 
s > 0, we can find a solution order by order. On the 
other hand if the determinant vanishes for some s, the 
equations can lead to inconsistencies. 

One remark has to be made about the order of mag- 
nitude of the coefficients. In the general case with 
four fields, we have found constraints on the coefficients 
An, Bn, Dn and from the observations. Since the dif- 
ferent fields are not related to the same observation, the 
constraints are different for each field. We have there- 
fore found that £)„ and En can be much larger than An 
and Bn for small n, see table (|ll|. For large n the situa- 
tion is reversed. Nevertheless, from the previous analysis 
of the equations, we see that this situation is not satis- 
fying. Indeed, if two coefficients are much larger than 
the two others at order s, it means that we can neglect 
the two small coefficients. Therefore, we introduce four 
new equations at order s but only two new coefficients 
and the determinant (|D23[) vanishes. Hence we can not 
insure that a solution exists. Therefore, even if the ob- 
servations allow less stringent constraints on the order 
of magnitude of some coefficients, the equations of mo- 
tion generically imply that all the coefficients can have at 
most the order of magnitude of the smallest one at each 
order. 



3. Application to our problem 

We need to transform the four equations of motion to 
apply the method described above. First, we multiply 



each equation by the correct power of e'^ and such 
that each power in the equation (jD3[) is positive or null. 

After this modification, equation ([7]) and (|13p have the 
correct form and the method can be applied directly. For 
n = 1/2, which is the case we consider in detail, equation 
([T^ reduces to 



e + v' e 



(D24) 



We can expand the two terms / = + ^ + 

and g = \fK, in power of a;, using the method described 
above. Let's call p the lowest order of the development 
of /, containing only the coefficients ai, 6i, d\ and ei and 
a the lowest order of g. 

At lowest order the equation becomes 



We have three possibilities: 



(D25) 



j^f^ ~ and g'^"^ ^ 0. The following order is then: 
j(p+i)^(<T) ^ Since g^"^' ^ 0, it implies = 
0. The same development can be made at each 
order, and therefore we find / = 0. 



• 7^ and g'^"^ — 0. As previously this implies 
g = 0. This case is the trivial case K, = and is 
therefore not interesting. 

• /'''^ .9'"' = 0. 

In this case, the following order is /(p+i) — 0, 

which has exactly the same form as equation (jD25p 
and implies therefore either / = or 5 = 0. 

Since we want /C 7^ 0, equation (fT^ becomes 



/ = 



C + C' 



= 



(D26) 



at each order. We recover one of the Schwarzschild 
equation for which the method can be applied easily. 
Note that for the three equations d?]), ^ and (|D26|) . 
the relations between the usual coefficients of Post- 
Newtonian parametrization are not modified by the addi- 
tional coefficients up to an order 10~^ which is the order 
of magnitude of the largest additional coefficient. 

Equation (jlip is slightly different. Indeed, the left- 
hand side is proportional to /C"+^, that means pro- 

/ \2(«+l) 

portional to ( I , whereas the right-hand side 
is proportional to {^^^i^ ■ So the left-hand side is 



\ Mr, J 



2n 

~ (10^'^)^" times larger than the right-hand 

side. Therefore, this equation can modify the relation be- 
tween the usual coefficients. Indeed, we will now mix the 
usual coefficients coming from the right-hand side, with 
the additional coefficients of the left-hand side which are 

2n 
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We consider in the following the case n = 1/2. The 
lowest order of the left-hand side is 8 whereas the lowest 
order of the right-hand side is 5. At order 5 equation 
(fTTj) will then have the following form 



1 



-F{Ai,Bi,Di,Ei,ai,bi,di,ei) = G{ai,bi, a, di) 



Mrsa 

(D27) 

where F is proportional to Ai,Bi, Di and Ei . This 
equation gives a relation between the usual and the ad- 
ditional coefficients. 

Orders 6 and 7 have the same form, except that they 
contain also the coefficients 02,62... and Dq. Order 8 is 
different. Indeed, at this order the left-hand side con- 
tains also ai,6i,di and ei. Since they are multiplied by 
10^^, we can neglect the terms coming from the right- 
hand side. The same occurs for all the following orders. 
Therefore, for orders 8 and larger, equation pTjl becomes 

+ 2^^ + /C"+i = . (D28) 
The same argument as in Eq. (ID24P implies that 



+ 2^ + V I = 
r 4 ' 



(D29) 



This is the second Schwarzschild equation, but it is valid 
only for orders larger than 7 in the development in power 
of X. For smaller orders we have to take into account the 
right-hand side. 

The same situation occurs for powers smaller than 5. 
The left-hand side can be neglected with respect to the 
right-hand side, and therefore we can apply the method 
to equation (|D29p . 

To summarize, the method can be directly applied to 
equations ^ and ((D26)) which replace Then 
we have to solve orders x'^,x^ and x"^ of equation pT|) . 
Our method can then be applied to equation (jD29[) for 
orders larger than x^ and smaller than x^. 



4. Solutions 

Let us divide each equation by the correct power of x 
such that the lowest order , containing only ai, 61, di and 
ei is X for each equation. Concerning equation (|lip . this 
means that we have to divide by x^. 

At order x for the equations ([7]), (|13p and (|D26p we 
find: 



Order x of equation pT|) will be treated separately 
since it contains also Ai,Bi,Di and Ei. 

We impose oi = — 1 to recover Newton's theory and 
therefore we find 



61-1, 
ei = 1/2. 



(D31) 



From the last equality of equation (|D30p we have two 
possibilities: either di = or C3 = — ci. 



a. di = 



At order x equation ^3)1 implies 



Hence, either C3 



(ci + C3)d2 = . (D32) 
-ci or d2 ~ 0. At order x^ we have 



C2 - (ci + C3) ds = 0. (D33) 

Again we have two possibilities: either C2 = ci 4- C3, or 
ds = 0. The situation is the same at each order. Indeed if 
we assume that we have chosen di = d2 — ... — dg-i — 
from the order x to x''~^, order x" gives 

(ci+C2+C3)s^-3(ci+C2-KC3)s-2(ci-C2+C3) 4=0. 

(D34) 

This means that the only possibility to have at least 
one of the dg 7^ 0, is to satisfy one of the relation 

(Ci +C2+ 03)5^ - 3(ci -I- C2 + 03)5 - 2(ci - C2+C3) ^0 , 

(D35) 

for some positive integer s. Therefore, we have to con- 
sider two situations: either dg = + 0{Do, Eq) for all s, 
or one of the relation (jD35[) is satisfied. 

In the first case, the usual parameters ds are equal to 
zero up to the order of magnitude of the additional coef- 
ficients which of course modify equation (|D34[) . We can 
calculate the positive order x,x'^... including the addi- 
tional coefficients. And we have also the negative order 
x'^ ,x~^ , .... We can show that these sets of equations 
imply either 



ds ^ Dq , Vs and Dq ^ Di ^ 1)2 •■ 
or 



(D36) 









di ^ Eod2 ^ E^ds ^ ... r 




bi = 


-ai , 




ds Eods+i ^ EQds+2 ^ ... ^ 


Eq 


ei = 




Do - Eodi r 


^ -Eq doo 


= 


(ci -1- cs)di . 


(D30) 


Ds ^ Eobs-i ^ ... r 





(D37) 
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Since -Doo ^ in order that the expansion ([2T|) con- 
verges in the Solar System, the first case imphes (r) = 
0. 



In the second case, since 



and dno is finite, we 



also have A'^{r) = 0. 

Hence A* (r) becomes the only degree of freedom of the 
vector field, which is then completely determined by the 
constraint d?]). 

In the following we will study in details the second 
situation where one of the ds at least is different from 
zero. We consider the simplest case where di ^ and 
therefore C3 = — ci. 



C3 



-ci 



From order of the four equations ©, ((T3|), p26p 
and (|D29|) . and using eq. (jD3ip wc find 



0-2 = 
62 = 

= 2c2df - 4c2(ii + ci - ■ 



4 

2 ' 



1 



(D38) 



We see that bi — 1 and 02 = 1/2 are in complete 
agreement with the observations. The additional coef- 
ficients imply only a contribution of order 10~^, which 
we have safely neglected since they are well beyond the 
precision of the measurements. The last equation allows 
to calculate di as a function of ci and C2. 

We have to determine the equations for orders and 
x"^ which come from the mixed terms of equation (|11|1 . 
We use the solutions for the previous coefficients. 

From order x^, we then find 



03 

63 
d3 



3 

^16' 
1 

16 ' 

— + -dl + did2 , 

( - 8d2C2 - 4cidi - 24cid? + 52dlc2 
0C2 \ 



+4:Cid2 + C2C?i -I- 32(i2C2(ii 
From the order we obtain 



(D39) 



1 



64 



,,,Jc2 + 80c2dl 
-l-816c2dt - GAcidl 



6Adid2C2 + 6AdiCid2 
384ci4 + 512d2C2d? 



-64d^C2 



192c2 



104cid2C2 + 70cidiC2 - 208cidiC2 



-16c?d2 - 16cldi - 32cid2C2rf? - 96c^d? 
-3264cid^c2 + 132d2cl - 652d?c^ - Tdict 



384d2C^di + 3840c^d:jd2 + TlOAcid 



2j4, 



„2 ,5 



-384c2did2 



(D40) 



We also need to solve for the negative powers. Using 
the hierarchy between the additional coefhcients \Ai \ ^ 
1^2]... we find from equations ([7]), and (jD26[) at order 



^0 = 
ii = 
A2 = 



and x that 



3^2 

16diC2 



24c2di - Scidl + Uc2dl - \\c2 - 2c\ 



B2 , Bi — 
—B2 , E2 



5B3 
4 ' 

2 



D2 



B_2 

2 

B2C1 



3^2 

4(ii 



8c2di 



(D41) 



We remark that Bi is of the order of B^, therefore it 
can be neglected with respect to Ai , Di and Ei , but also 
with respect to the order 2. B2 remains undetermined. 
Bj, will be determined by the lower order x~^ , but we are 
not interested in its value here. 

We now consider equation (|lip at the order x, x^ and 
x~^ , which contains both the coefficients of negative and 
positive powers. We solve this system of three equa- 
tions plus the equation (jD30p with the help of maple. 
It contains the five variables di, ^2, ci, C2 and B2 and the 
two parameters ai and M . We find a set of solutions, 
from which we only consider those with ci and C2 nega- 
tive in order to ensure a positive-definite Hamiltonian for 
perturbations and non superluminal propagation of spin- 
degrees of freedom in the approximately-Minkowski 
regime of the theory [28| 



ci = 
di - 

Cl 



-33.11 / B. 



0.16, 



\MrJ 



C2 



-30.75 / B2 n2 

a? \Mr, 



dn ^ -0.62. 



(D42) 



0.48/ B2 \2 



\Mr, 



C2 



6.36/ B2 \2 



'-( 



-0.10, 



0.28, 



(D43) 
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Cl 



di = OM, 



-140.29 / B2 \2 



\Mr, 



C2 



61.93/ B2 ^2 



a 



do = -7.89. 



(D45) 



Each solution implies a strong constraint on the pa- 
rameters ci,C2 and ai. Indeed B2 = -B2 (7^) '^^'^ 
at most 10~^2 in order not the be detectable in the So- 
lar System, and therefore (y-^^^ — 10~^*. This means 
that either ci and C2 or ai has to be very small if ai ^ 0. 

Finally, we have to calculate the determinant of the 
system to determine if a solution exists at each order. 
Using the method described above we find from the pos- 
itive powers in x (which correspond to negative powers 
in r){s > 1) 



/ as ^ 




( \ 


hs 






ds 













where 



(D46) 



/ -1 







5 = 



— s 

dis{s—l) 
2 

r(s) 



-2 \ 



ci(s-l)(s-2) rfis(s-l) 
4 2 



s(s-l) 

-S 



(D47) 

and Hi"' are functions of the previous coefficients 
oi, ...a^-i, 5i, di, .., 4-1, ei, .., e^-i. Therefore, 
at order s a unique solution (os, 6^, dg, eg) exists if 

det5 = -ys2(s-l)2(s-2)7^0 . (D48) 



Since ci 7^ 0, the condition is satisfied for s > 2. The 
order s = 2 for which the determinant vanishes has been 
solved explicitly above and leads to no inconsistency. The 
solution is not unique, since we have found four different 
values for ^2- For each order (s > 2) the determinant 
is non zero and therefore we can determine the unique 
solution as a function of the previous orders. 

The same structure repeats for positive powers of r 
is>2) 



S 



S = 



(aA 






Bs 






Ds 






\Es) 






1 












-s{s + l) 



where (D49) 



cis{s-\-l)di 
2 



is) 

and are functions of oi, 61, di and ei, and 



C2(14-s)d2- ci^(s+l)di 
ci(s-l)(s-2) 2 

(D50) 



,B 



s-l, 



of the previous coefhcients Ai,...,As-i, Bi, 

bi,...,bs-i,Ei,-,Es-i, and Q^"' < h['\ A unique 
solution exists at order s if 



detS" = -2s2(s + i)(^c2(14-s)d?-^(s-l)(s-2)) ^ . 

(D51) 

The orders s — 0,1 and 2 have been calculated above. 
For each of the solutions p42|) to (|D45|) we find that 
det 5 ^ V s. 
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